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Abstract 

The Bethe-Salpeter equation provides the most widely used technique 
to extract bound states and resonances in a relativistic Quantum Field 
Theory. Nevertheless a thorough discussion how to identify its solutions 
with physical states is still missing. The occurrence of complex eigenval- 
ues of the homogeneous Bethe-Salpeter equation complicates this issue 
further. Using a perturbative expansion in the mass difference of the 
constituents we demonstrate for scalar fields bound by a scalar exchange 
that the underlying mechanism which results in complex eigenvalues is the 
crossing of a normal (or abnormal) with an abnormal state. Based on an 
investigation of the renormalization of one-particle properties we argue 
that these crossings happen beyond the applicability region of the lad- 
der Bethe-Salpeter equation. The implications for a fermion-antifermion 
bound state in QED are discussed, and a consistent interpretation of the 
bound state spectrum of QED is proposed. 
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1 Introduction 



In a relativistic Quantum Field Theory bound states and resonances are identified 
through the occurrence of pole or cut like singularities in the Green's functions of the 
theory. Thus it is natural to study the four-point Green's function when searching for 
two-particle bound states. Indeed, the corresponding equation has been proposed by 
Bethe and Salpeter M and subsequently been proven by Gell-Mann and Low |2| and 
Schwinger || as early as 1951. It can be interpreted as one of the Dyson-Schwinger 
0. Q equations of Quantum Field Theory and it is, thus, an inhomogeneous integral 
equation. Assuming the existence of a bound state signaled by a pole in the four-point 
function the homogeneous Bethe-Salpeter equation may be derived (see below). It 
allows to determine the bound state masses and covariant wave functions. Usually one 
employs the so-called ladder approximation which renders the homogeneous Bethe- 
Salpeter equation in the form of an eigenvalue problem. Hereby the eigenvalue is the 
square of the coupling constant, and the bound state mass has to be tuned such that 
the Bethe-Salpeter eigenvalue equals the given value of this constant. The covariant 
wave functions are then determined as the eigenfunctions of the system. 

The only analytically solvable example of a Bethe-Salpeter equation is the one for 
two (massive) scalar particles bound by the ladder approximation to the exchange of 
a massless scalar field || |6|, |?J . Despite its relative simplicity as compared to realistic 
systems this model, the Wick-Cutkosky model, displays already the advantages (e.g. 
full covariance) as well as the shortcomings (e.g. the existence of abnormal states) 
inherent to almost all Bethe-Salpeter based approaches used until today, for a recent 
very accurate numerical solution of the ladder Bethe-Salpeter equation for scalars 
see e.g. || and references therein. 

Common to the analytical and the numerical solutions is the existence of ab- 
normal states which have led to controversial discussions regarding their physical 
interpretation P, |5|, [7|, [Sj, [TQfl. Even worse, for the case of constituents with unequal 
masses some eigenvalues of the homogeneous Bethe-Salpeter equation become com- 
plex [IT], |13|]. Clearly, such a behavior is unexpected and has to be understood. It 



is usually attributed to the use of the ladder approximation which destroys crossing 
symmetry from the very beginning. However, we will see in the following that already 
the use of bare propagators is problematic. This will also not be cured by incorpo- 
rating one-particle self-energies and thereby generalising the ladder approximation. 
Instead, it will be shown in this paper that the occurrence of complex eigenvalues 
also happens with self-energies taken into account. It seems that one has to start 
on a more fundamental level: One seriously has to raise the question which values 
of the renormalized coupling constant are possible when considering a renormalized 
Quantum Field Theory. 
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1.1 Derivation of the homogeneous Bethe— Salpeter equation 

In order to make this paper as self-contained as possible, and to make our argumen- 
tation accessible to non-expert readers, we will supply some important derivations 
and facts concerning the Bethe-Salpeter approach to relativistic bound states in the 
remainder of this introduction. Readers familiar with the Bethe-Salpeter equation 
probably will probably prefer to jump to Sec. 2 immediately. 

We assume to deal with three types of scalar fields. The two constituents are 
supposed to have masses mi and m 2 and self-energies Ei and E 2 . The four-point 
function (x±, x 2 ,yi, y 2 ) describing the scattering of these two constituents fulfills 
the inhomogeneous Bethe-Salpeter equation 

[D X2 +m\- E 2 ] [D X1 +m\- Ex] G (4) {x u x 2 , Vl , y 2 ) = (1) 

5 (4) (xx - yi) 5 (4) (x 2 - y 2 ) + 6® (x x - y 2 ) 5 (4) (yx - x 2 ) 

+ J d A z x d A z 2 K{xx y yx,z 1 ,z 2 )G { ^ (z t , z 2 ,y t ,y 2 ) . 

The kernel K is defined as the sum of all amputated two-particle irreducible contri- 
butions. The Feynman diagrams of the first few terms are depicted in fig. [1|. 

For translationally invariant systems it is, of course, advantageous to transform 
this equation to momentum space. Note that the introduction of the relative coordi- 
nate x = Xi — x 2 allows an arbitrary parameter a G [0, 1] in defining the coordinate 
X = ax\ + (1 — a)x 2 which results in the corresponding momenta 

Pi = aP + p and p 2 = (1 — a)P — p (2) 

for the constituents in terms of the total and relative momenta, P and p, respectively. 
Fourier transforming eq. (JJ) leads to 

f-pL[D (p, p', P) + K (p, p', P)} G (4) (p', p", P) = 5^ (p - p") , (3) 

where 

D(p,p',P) := {2k) a 5^{p-p>) (G^y'ipx) (G^y\-p 2 ) (4) 

is defined in terms of the inverse two-point Green's functions of the constituents. 

The crucial step from the inhomogeneous to the homogeneous Bethe-Salpeter 
equation consists in assuming a bound state reflecting itself in a pole in the four- 
point Green's function for an on-shell momentum P os , with P^ s = M 2 

G (4) (P,P', Pos) = ~\*!?f; S ^ {p! : P i + reg. terms, u := V^T^P (5) 
(2ix) 2u (_r u - uj + ie) 
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Figure 1: Feynman diagrams of the first few terms in a perturbative expansion of 
the kernel K. Solid lines represent propagators of the constituents, dashed lines the 
propagator of the exchange particle. 




Figure 2: Pictorial representation of the homogeneous Bethe-Salpeter equation in 
ladder approximation. 



where we introduced the definition of the Bethe-Salpeter amplitudes 
X (x 1 ,x 2 ,P) := (0|T{$(zi)$(z 2 )}|P) 

X (x 1 ,x 2 ,P) := (Q\T{&( Xl )&{x 2 )}\P) (6) 
together with their Fourier transforms 

x(x u x 2 ,P) =: e~* px [ -^e-^ x (p,P) 

X(X!,X 2 ,P) =: e+* px [ 7^je-*»x (p, P) ■ (7) 

J (2vr) 

Hereby |0) denotes the ground state (vacuum) and \P) the bound state. Very close 
to the pole the regular terms can be safely neglected and the dependence of the four- 
point function on the relative momenta p and p' can be separated. Expanding 
and [D — K\ in the inhomogeneous Bethe-Salpeter equation in powers of (P° — uo) 
yields the homogeneous Bethe-Salpeter equation and the normalisation condition for 
the amplitude. The order (P° — lu) 1 provides 

/ 44 i D (P>P'> P os) + K {p,p\ P os )\ x (p, Pes) = 0, (8) 
J (2tt) 



whereas to O [{P Q — co>)°j one obtains 

/ 7^7^ tr (* <P> P ^ 7^-o ( D p ) + K M p )) X (p', Pas)) 

J (Ztt) (An) V or p0=u) J 

= 2iu. (9) 

This ensures the residue to be equal to 1 at the bound state pole. 

The homogeneous Bethe-Salpeter equation (|3|) is a linear integral equation for the 
amplitude \ whose overall normalisation is fixed by (||). Approximating the kernel 
by the one-boson-exchange depicted in the first diagram of fig. eq. flD can be cast 
into an eigenvalue problem for the coupling constant by defining the vertex function 

T( P ,py. 

x(p,P)=:G 1 (p 1 )G 2 (p 2 )T(p,P). (10) 
The homogeneous Bethe-Salpeter equation (§) then reads 

r (p, p os ) = - [ p^K ( P ,p', p os ) g x (pi) g 2 ( P ' 2 ) r (j/, p os ) . (ii) 

In the ladder approximation the kernel K is set equal to 

K(p,p')=i- ^ (12) 

(p — p) 1 — fi z 

which is just the propagator of the scalar exchange particle of mass [i multiplied with 
g 2 . On inspection one finds that (|TTD is an eigenvalue problem for g 2 if G\ and G 2 are 
the bare propagators of the constituents. The ladder approximation to the Bethe- 
Salpeter equation is pictorially represented in fig. |2|. If a parameter pair (g 2 , P° = M) 
exists the pole assumption is a posteriori justified and M is the bound state mass 
with x being the corresponding amplitude (wave function) as can be inferred from 
eq. (H) which reflects, of course, nothing else than the Lehmann representation of the 
four-point function. 



1.2 Abnormal solutions and relative time parity 

For a given bound state mass M the eigenvalue spectrum of the Bethe-Salpeter 
equation should be positive definite, i.e. the eigenvalues should be real and positive. 
We will see that this is not the case. Furthermore, one expects that for very small 
coupling constants the binding energy vanishes and the bound state mass becomes 
identical to the sum of the masses of the constituents. There are states, however, 
which possess vanishing binding energy for a finite coupling. They are therefore 
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Figure 3: The eigenvalues of the Bethe-Salpeter equation in the Wick-Cutksoky 
model with constituents of equal masses as a function of the ratio of bound state 
mass to the sum of masses of the constituents,^ = M/2m . 




called abnormal states. In the Wick-Cutksoky model [^], [5| with constituents of 
equal masses mi = m 2 = m these abnormal states are easily identified: They only 
exist for A := g 2 /lQ7T 2 m 2 > A c = 1/4. If the binding energy becomes very small, 
i.e. t] := M/2m — > 1, the corresponding coupling constant A vanishes for the normal 
solutions, i.e. A — > 0, whereas A — > A c = 1/4 for the abnormal states, see fig. § and 
its closeup fig. £|. 

On the other hand, in the opposite limit of a massless bound state, i.e. for 77 — > 0, 
the eigenvalues are given by 

A = (n + K)(n + K + l) =N(N+1). (13) 

The value of k equals the number of nodes of the eigenfunctions when plotted as a 
function of the relative time Xq. For all normal solutions, including the ground state, 
one has k = 0. This especially means that all normal solutions are even under the 
reversal of relative time Xq — > —Xq. On the other hand, k — 1,2, . . . leads to abnormal 
solutions, i.e. these solutions may be even or odd with respect xq —>■ —xq: 

x(-x Q ,x) = (-l) K x(^o,x) . (14) 

One may interpret the abnormal solutions as excitations in relative time. They will 
obviously not appear in a purely non-relativistic treatment where the constituents 
are considered for equal times only. However, not all abnormal solutions necessarily 
vanish in a three-dimensional reduction of the Bethe-Salpeter equation fl4j| . On 



the contrary, the spectrum of a three-dimensionally reduced equation will contain 
remnants of these abnormal states. 

For constituents of unequal masses and a massless exchange particle the Bethe- 
Salpeter amplitudes are asymmetric, i.e. neither even nor odd, under the transfor- 
mation x — > — Xq. However, the eigenvalues of the abnormal states start all at some 
common A c , see also |7j and references therein. 

For equally massive constituents and a massive exchange particle the amplitudes 
do have a definite "parity" under the inversion of relative time. However, there is no 
longer a common A c at which the eigenvalues of the abnormal solutions start. 

The controversial discussion of the abnormal states already started with the work 
of Wick |J and Cutkosky |]] , for some of the related investigations see also || [T5[ |7|] . 



Fixing the spatial coordinates of the constituents leaving only the relative time as 
degree of freedom one is able to show that the abnormal solutions are all unphysical 
for this "static" model flHj. However, it is not clear whether abnormal solutions can 



in general be considered as unphysical [|16|]. As already mentioned, some of these 



abnormal solutions survive most three-dimensional reductions [Q. Based on the 
analysis of several reductions of the Bethe-Salpeter equation the main conclusion of 
ref. mj is that abnormal solutions are very probably spurious consequences of the 



ladder approximation supporting hereby an old conjecture by Wick 
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Figure 5: The lowest lying three eigenvalues of the ladder Bethe-Salpeter equation 
for 5 = m 1 /m2 = 4 and ji = m 2 . 



M,fi 


symmetry group 


Bethe-Salpeter amplitude 


M = 0,fi = 


0(5) 


X oc Z nMm (fi 5 ) 


M = 0,/i ^ 


0(4) 


X oc Z fe/m (Q 4 ) 


M/0,|i = 


0(4) 


X 4m (^4) 


M ^ 0,/i ^ 


0(3) 


X oc y /m (Q 3 ) 


/! — > OO 


0(4) 


x oc z klm (n 4 ) 



Table 1: Summary of the symmetries of the scalar Bethe-Salpeter equation in ladder 
approximation, ji denotes the mass of the exchange particle and M is the Mass of 
the bound state. The functions Z (or Y) denote the spherical harmonics for the 
corresponding n-sphere Q n . 
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1.3 Complex eigenvalues 



In this subsection we will consider the general case for a Bethe-Salpeter equation 
containing only scalars: the masses of the constituents are assumed to be unequal, 
the exchange particle is chosen to be massive. The numerical method is detailed in 
appendix |A| and we will only give results here. For comparable sets of parameters our 
results are equal to those given in refs. [T3| . In fig. ^ the lowest three eigenvalues 
for 5 = mxjm^ = 4 and \i = m 2 are shown as a function of rj = M/(mi + m<i). 
The eigenvalue of the ground state is real for all physically allowed values rj e [0, 1] 
and vanishes for 77 — > 1. The ground state is thus a normal state. The situation 
differs drastically for the higher-lying two states. At r\ — 1 the eigenvalues are non- 
vanishing, real and positive. As the binding energy increases these two levels become 
degenerate at 77 « 0.81. In the interval given approximately by 0.25 < 77 < 0.81 the 
eigenvalues are complex with a degenerate real part and imaginary parts of opposite 
sign. As one can see in fig. |5| from the projection of the curves not all eigenvalues are 
monotonically decreasing as a function of 77. Clearly, such a behaviour is unphysical 
because a decreasing coupling constant should result in less binding. 

The appearance of complex eigenvalues is not restricted to a special choice of 
parameters. In fig. ^| the absolute value of the imaginary part of the eigenvalues 
for the first two excited states is shown for fi = 1712 and 1 < 8 < 15. One clearly 
sees that with the increase in the mass of one constituent the interval in 77 for which 
the eigenvalues are non-real becomes smaller, however, the maximum value of the 
imaginary part even increases. For a given mass ratio 5 complex eigenvalues exist for 
masses of the exchange particle up to some 7i max , i.e. for < ji < /U max - Considering 
the first two excited states we estimate this 7i max to be 7i max ~ 0.85m! + 1.66m 2 . 
Increasing the mass of the exchange particle the higher lying eigenvalues tend to 
become real again. This can be understood from the fact that for an infinitely heavy 
exchange particle the Bethe-Salpeter equation assumes an 0(4) symmetric form as 
in the case of an massless exchange particle, see also table [I] which summarises the 
symmetries of the scalar ladder Bethe-Salpeter equation. 



2 Inconsistency of the abnormal states 

In this section we relate the appearance of complex eigenvalues of the Bethe-Salpeter- 
(BS) equation to the presence of the abnormal solutions. To this end we verify a 



conjecture of Kaufmann [] I 



In section |1.1| we recalled the steps leading from the inhomogeneous BS equa- 
tion (P to the homogeneous BS equation (§). We now Wick- rotate the contour of 
integration in (H) according to po — * ip±, P — * P an d thus obtain the homogeneous 
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Figure 6: The absolute value of the imaginary part of the eigenvalues for the first 
two excited states. 

BS equation for euclidean momenta. A subsequent Fourier-transformation yields the 
homogeneous BS equation in euclidean coordinate-space: 



The operator C is the product of the inverse propagators of the constituents and \ is 
the BS amplitude describing a bound state of two scalar particles with masses mi and 
m 2 , respectively. The interaction of the two constituents is encoded in V which in 
the ladder-approximation is simply the Fourier-transform of the bare propagator of 
the exchange-particle. We decomposed the operator C into two parts: Co, containing 
the terms that are present for mi = m 2 and V containing the additional terms that 
appear only for m\ ^ m 2 . The eigenvalue is denoted by A instead of A for later 
convenience. The explicit expressions are given by: 



C X =[C + V}x = AV X . 



(15) 




(16) 



with p^ = id^/m 2 



V 



+ 



(ip 4 ) [4 V A [p 2 -(l 
A 2 [2 (1 - „ 2 ) [j? - 



(1 - V 2 )] ~ Vp2] + A 4 (1 - r?) 2 



(17) 



V 




(18) 



A 



167r 2 m2 



(19) 
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The parameter rj = M/(m% + m 2 ) measures the mass of the bound state, A = 
(mi — m 2 )/2m 2 is proportional to the mass difference, and \i is the mass of the 
exchange-particle. 

The explicit expression (|17|) shows that V vanishes for a vanishing mass difference 
A. Thus for small A the decomposition C = Co + V can be understood as dividing 
£ into the dominant part Co and into the perturbation V. In the following we will 
formalise this point of view and investigate the consequences thereof. 

Suppose Xi an d X2 are two solutions of the homogeneous BS equation for con- 
stituents of equal mass with eigenvalues Ai and A 2 respectively: 

C0X1 = AiVxi 

£0X2 =A 2 Vx 2 (20) 

Because of the restriction mi = m 2 the BS amplitudes xi an d X2 have a definite 
xo-parity. As we will see the interesting case is to take Xi an d X2 to be of opposite 
a^-parity, and this is assumed for the following considerations. Since we are interested 
in the region where two real solutions become a pair of complex conjugated solutions 
we take the ansatz x — a Xi + ^X2 where x is the BS amplitude for different masses 
mi 7^ m 2 and Xi, X2 are the amplitudes for mi = m 2 as defined in (ppp. Substituting 
this ansatz for x m t° (0) an d projecting onto xi an d X2 one obtains a system of 
equations relating the eigenvalue A to the eigenvalues Ai, A 2 : 

aAiVn + aVu + bVn = AaVn 

6A 2 V 22 + aV 2 i + bV 2 2 = A6V 22 (21) 

where 

X lk = ( X i, X Xk ) ■= I d A x x* (x) X Xk (x) (22) 



for X = V, V and i, k = 1, 2. In deriving (^TJ) we used V12 = V21 = which is due to 
the different x -parity of Xi an d X2- 

Assuming xi an d X2 to be normalised according to Vn = V 22 = 1 one finds the 
two eigenvalues Ai 2 of the BS equation (O), 



Ai )2 = \ (Ai + A 2 + Vn + V 22 ) ± ^ (Ax - A 2 + P n - V 22 f + P 12 P 2 i . (23) 

This gives A, i.e. the eigenvalue for A 7^ 0, in terms of the matrix-elements Vik and 
Ai, A 2 which are the eigenvalues for A = 0. Using the explicit expression QP7f) one may 
calculate the matrix elements Vik thereby employing the solutions xi an d X2 obtained 
by numerically solving the BS equation for A = 0. On inspection one finds that V is 
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Figure 7: The real part of two colliding solutions, according to the numerical solution 
(solid line) and its semi-analytic estimate ( [23] ) (dashed line), c.f. also the upper left 
corner in fig. [j| 

antihermitian leading to Vn = —V^i => V\%Pi\ = — \Pvi\ ■ According to (^) this is 
a necessary condition for the appearance of complex eigenvalues. This is in contrast 
to the hermitian perturbations that are considered usually. Whereas a hermitian 
perturbation will enforce a repulsion between the solutions (avoided crossing) an 
antihermitian perturbation will result in an attraction. 

The real part of Ai,2 as calculated according to fl2"3"| ) together with the result of 
the numerical calculation is shown in fig. [j]. The dashed line represents the results 
according to the coupling mechanism that has been proposed in this section whereas 
the solid line represents the results of the numerical solution of the BS equation with 
5 = mi/m 2 = 4 and (1 = 1712, c.f. the upper left corner in fig. ||. The two curves 
are in reasonable quantitative agreement. The simple ansatz x — a Xi + ^X2 and the 
coupling due to V reproduce the overall appearance of the collision astonishingly well. 
The imaginary parts of the solutions A (77) appear after the collision point. They agree 
with those of the numerical solution to the same level as the real parts. Therefore we 
conclude: the coupling of solutions of opposite xo-parity is the mechanism leading to 
complex eigenvalues. Hereby, the presence of solutions of positive and negative xq- 
parity is a prerequisite for the appearance of complex eigenvalues A. Choosing both 
states to be of the same x -parity will lead to an avoided crossing. Thus, at least 
one abnormal state is necessary for the occurence of a pair of complex eigenvalues. 

In this context it is interesting to note that Naito and Nakanishi ([HJ derived a 
normalization condition which they claimed to remove the solutions with complex 
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eigenvalues. However, a substantial loophole has been spotted in their derivation 



and Fukui and Seto |T3[ showed numerically that the normalization condition can be 
satisfied for solutions with real as well as for solutions with complex eigenvalues. 

In the following we will show that the solutions with negative xo-parity lead to 
inconsistencies on a even more fundamental level. We will derive an equation that 
has to be satisfied by all solutions of the BS equation in ladder approximation which, 
however, is only satisfied by the solutions with positive x -parity, see also refs. j[9, . 

The equation may be derived starting with the inhomogeneous BS equation 

D (p, P) G (p, q, P) = 5 A (j>-q)+ J d A p K (p, p', P) G (p', q, P) (24) 

where D is the product of inverse propagators 

D(p,P):=G 1 (p 1 )G 1 (p 2 ), Pl :=aP + P , p 2 := (1 - a)P - p } (25) 

and is slightly different from D defined in ([!]). The following considerations are 
restricted to the case m\ = m2 which is reflected in the definition (p5|) . G is the 4-point 
function. It depends on the total 4-momentum P and the relative 4-momenta p and 
q. The kernel K is defined as the sum over all 2-particle irreducible diagrams. Taking 
the derivative with respect to A = g 2 / lQir 2 m\ one obtains after some reordering 



d A pdyG{p",p,P) 



—D {p, P) 5 4 {p' -p)-q^K (p, p', P) 



G (p', q, P) 



A G (p",q,P) . (26) 



No approximation has been used so far. To proceed we restrict the considerations 
to the case of bare propagators and to the use of the ladder approximation for the 
kernel K. This yields 

= o, 

^K(p,p',P) = ±K(p,p',P). (27) 

Substituting the representation of the 4-point function as given in eq. (|5|) and com- 
paring the pole contributions one obtains 

r dM 2 
i J d 4 p d 4 p' x (p) K (p, p\ P) x ip') = A^ (28) 

where x is the BS amplitude. Equation ( p8|) should be satisfied by all solutions x °f 
the BS equation if bare propagators and the ladder approximation for the kernel K 
are used. 
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One can make use of translation invariance to factor out the center of mass motion: 



x = e- ip - x <f> (x, P) . (29) 

The coordinates are hereby defined according to 

x := X\ — x 2 , X := a.X\ + (1 — a)x 2 (30) 

and the Xo-parity is reflected in the relation (— x°, x) = e (0) (x°, x) with e (0) = 
±1. In the next step one decomposes the amplitudes into positive and negative 
frequency parts, see ||. After Wick rotation, 

V := (Po, P) ->• (*P4, P) =: P, (P 2 > 0), (31) 
one obtains the relation, see appendix [FJ for its derivation, 

0(p,P)=e(0)[0(p,P)r (32) 

for the Euclidean amplitudes. Applying the Wick rotation (|3T| ) to and substi- 
tuting the relation (|32| ) one arrives at 

- e (0) J d'pd'p' [0 (p, P)]* K' (p,p>, P) (p', P) = (33) 

with positive K' (p,p', P) = (p — p') 2 . As was discussed in the introduction a 
decreasing coupling constant should result in less binding, i.e. dM 2 /d\ should be 
negative. Eq. (|33|) implies that this is the case only for e (0) = +1, i.e. only the 



states with a positive Xo-parity can satisfy both, eq. (|33D and the physical condidtion 
dM 2 /d\ < 0. 

On the one hand, we have obtained eq. (|28|) that has been derived from the exact 
BS equation using bare propagators and the ladder approximation for the kernel. On 
the other hand, we found that the solutions of negative Xo-parity (being solutions of 
the BS equation in the same approximation) cannot satisfy this equation and behave 
simultaneously reasonable under the change of the coupling constant. As a matter 
of fact, we found in our numerical results examples for both types of behaviour, i.e. 
either a violation of eq. ( |33"D or a positive slope when ploting M 2 vs. A. As a subset 
of the abnormal solutions the solutions of negative Xo-parity are present only above 
a certain nonzero coupling-constant. To reconcile these conflicting results one has 
to question the combined use of free propagators and the ladder approximation for 
coupling constants that are above this critical value. 
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3 Renormalization of one— particle properties 



In this section we first present a consistent approximation scheme for the propaga- 
tors and the kernel of the BS equation. The resulting system of Dyson-Schwinger 
equations is solved yielding a domain of validity for the ladder approximation. The 
equations are derived using the formalism of Cornwall, Jackiw and Tomboulis (CJT) 
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The starting point is the generating functional 



Z[J,K] := J JJX^exp 



S[<S>] + $ • J + -$ • K ■ $ 



(34) 



where &i(x) are scalar fields, S[<&\ is the classical action and the abbreviations are 
defined by 



$ • J := [ d * x $ i ( x ) Ji ( x ) ' 



(35) 



The generating functional for connected Greens functions is defined as usual by 
Z [J, K] =: exp (iW [J, K)). Using 



6W [J, K] 

5J { (x) 
5W [J, K] 
5K ik (x,y) 



- {fc (x) <f) k (y) + G ik (x, y)) 



(36) 



and performing a double Legendre transformation one obtains the effective action T 
as a functional of d> and G, 



r [0, G] := W [ J, K] - ■ J - -0 ■ K • - -G ■ K. 
Here the notation defined in (|35|) has been used. In ref. []21| the expression 
T [0, G] = S [0] + Utt Ln (G- 1 ) + ~iTr (V^G) + T 2 [0, G] 



has been derived where T> 1 is given by 
iV ik ( x , V\ 4>) ■= iDr 1 (x - y) 8 ik + 



6 2 S int [0] 
5<pi (x) 5(p k (y) 



(37) 



(3* 



(39) 
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and D is the bare propagator. The nontrivial part in the effective action fl3"B|) is T 2 . 
It is defined as a sum of vacuum loops where the vertices are given by the interaction 
part of S[$ + 0] and the propagators are equal to Go-- 

The BS equation for the amplitude ip in its most general form is given by 

dV dV r£, fc , (*, y; 2/') (*', y') = (40) 

where V iki , k , (x, y; x , y) is the 2nd functional derivative of T with respect to G. The 
propagators G are the solution of the corresponding Dyson-Schwinger equations 

W -0. (41) 



5G ab (x, y) 



In the following we will consider a theory of three scalar fields $i,$2 and $3 
describing massive particles and interacting according to the Lagrangian 

£ = E 5 - E + 4 + ^2) $3. (42) 

i i v 

We approximate r 2 by the first nontrivial term of the loop-expansion, i.e. we use 

r ,-e 

= \g 2 j d'x'dS' (j2Gl(x',y')G 33 (x',y') (43) 

\i,k=l 
2 

+ E G ^ a;, ^ , )G'33(^,y , ) 
i=l 

This approximation for T 2 results in the rainbow-ladder approximation for the BS 
equation (50) and the Dyson-Schwinger equations (5T). Having evaluated the func- 



tional derivatives we Fourier transform the resulting system of equations to momen- 
tum-space thus obtaining the BS equation 

1-1 t„ \ n-i \ „/.i2 /„ d\ „2- f d q , . 12 



(Pi) Gj 1 (ft) ^ (p, P) =g<i * G 3 (p - g) ^ (g, P) (44) 

•/ ( 27T ) 

for the BS amplitude ?/> 12 describing a bound states of particles 1 and 2. The momenta 
are defined by 

P\-.— p — aP, p 2 :—p+ (1 — a) P. (45) 
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Figure 8: Graphical representation of the coupled system of Dyson-Schwinger equa- 
tions fl4"6|), (P7|). Thick and thin lines represent the dressed and the bare propagators, 
respectively. 



Evaluating explicitely the functional derivative fl4T|) we obtain the coupled system of 
Dyson-Schwinger equations: 

GJ 1 (p) = Di l (p) - 2g 2 i [JL«G a (p-q) G { (q) i = 1, 2 (46) 

7 (2tt) 

G 3 " 1 (p) = D 3 " 1 (p) / ?) ^ (g) (47) 

^ ( 27F ) i=l,2 

for the propagators Gj that are used in the BS equation. This system of equations is 
depicted diagrammatically in fig. (H). 

Up to now we have used the unrenormalized fields $j, masses rrii and coupling 
constant g. Introducing the corresponding renormalized quantities &i,rhi,g through 

\[Zi% = $ i5 Z m .m 2 = m 2 , Z g g = g (48) 

one is able to repeat the steps leading from fl41|) to fl46|) and fl47j). As we use 
bare vertices (rainbow approximation) we have to set Z g — 1 for consistency. This 
procedure yields the Dyson-Schwinger equations for the renormalized self-energies 

E« (p 2 ) = A t + Zflfi f^<G 3 (p-q) Gi (g) i = 1, 2 (49) 

7 (2tt) 

£ 3 (p 2 ) = A 3 + g 2 i [ -p- A z & (p - ?) Gi (g) (50) 

( 27r ) i=l,2 

where the abbreviation 

A, := (1 - Zi) p 2 + (Z t Z mi - 1) m 2 % = 1, 2, 3 (51) 
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Figure 9: The dependence of the renormalization constants Z^iji = m i) an d 
Z m2 {KjjVi H> — m 2) on A = g 2 /lQi[ 2 m\. The renormalized masses of the constituents 
were taken to be mi/1712 = 4 and the mass of the exchange particle was set equal to 

has been used. This coupled system of Dyson-Schwinger equations is finally defined 
by fixing the renormalization constants according to the on-shell renormalization 



Note that the Lagrangian ([^) defines a super-renormalizable model. Thus, the mass 
renormalization constants Z mi are logarithmically instead of quadratically divergent, 
the field renormalization constants Zi are finite. We take them nevertheless into 
account in order to mimic an exactly renormalizable theory as closely as possible. 

In order to obtain a numerical solution a Wick rotation to Euclidean space is per- 
formed. The resulting coupled system of integral equations is then solved iteratively. 
The details of the numerical method are given in appendix 0. 

Since we are interested in the domain of validity of the ladder approximation we 
focus on the physical constraints to be satisfied by the propagators. The quanta 
of the fields $ are assumed to describe physical particles. This implies that the 
associated renormalization constants Zi and Z mi are positive and smaller than 1, i.e. 



The dependence of the renormalization constants Zi and Z m . on the coupling 
parameter A = g 2 / IQ^m 2 , displays an interesting behaviour. Fig. |9] shows the 




_d_ 

dp 2 




0. 



(52) 



0<Z h Z mi < 1. 
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results for Z 2 and Z m , 2 which are the renormalization constants for the less massive 
field $2- The masses of the constituents were taken to be mi/m 2 = 4 and the 
mass of the exchange particle was set equal to m 2 . Z 2 and Z m2 start at one for 
small values of the coupling constant and decrease for greater values. Both, the 
mass renormalization Z m2 as well as the field renormalization Z2 become negative 
at some critical value of the coupling constant. These critical values are A ~ 0.062 
corresponding to g/m 2 ~ 3.13 for Z m2 and A ~ 0.086 corresponding to g/m 2 ~ 3.68 
for Z 2 . Since the propagator for particle 2 enters the BS equation we conclude that 
the domain of validity of the system of equations fl44|) , (f46|) , (|47D is limited to the 
range < g/m 2 < 3.13. A negative value for Z m2 especially implies that m 2 < 0, i.e. 
particle 2 becomes tachyonic. Note also that in ref. |22j the critical value A c ~ 0.063 
of the coupling constant was found using a variational approach. 

It remains to solve the BS equation thereby employing the dressed propagators 
which are the solution of fl49|) and fl50"|). Restricting the calculations to the allowed 
range of the coupling constant g one does find that the inclusion of the self-energies 
somewhat lowers the solutions g(M), however, the deviations are much smaller than 
1%. Leaving the domain of validity i.e. going beyond g/m 2 « 3.68 one does find a 
growing influence of the self-energies as could be expected. For gjm 2 ~^> 3.68 we also 
found complex eigenvalues A. 

The effect of the inclusion of the set of crossed ladders has been investigated by 
Nieuwenhuis and Tjon ||23|| . They used the Feynman-Schwinger representation which 
takes account of the ordinary as well as the crossed ladders. Within this framework 
they calculated the binding energy of the ground state and compared with the cor- 
responding results of the various bound-state equations. They concluded that the 
ladder Bethe-Salpeter equation substantially underestimates the binding energy of 
the ground state for large values of the coupling constant. This seems to support 
the results of this work. However, we caution the reader not to compare the re- 
sults directly. Nieuwenhuis and Tjon calculated the effect of the inclusion of the 
crossed ladders and used bare propagators within the Bethe-Salpeter equation as 
well as within the Feynman-Schwinger representation. This may be used to estimate 
the importance of the diagrams that have been neglected in the kernel of the ladder 
Bethe-Salpeter equation. The present work emphasizes the necessity to use consistent 
approximations for the kernel of the Bethe-Salpeter equation and for the self-energies 
of all particles. It has been shown that the renormalizability of the one-particle prop- 
erties restricts the domain of validity of the approximation that has been used for the 
effective action. Both approaches limit the applicability of the ladder approximation, 
but they do so using quite different constraining principles. 

To summarise: the physical condition < Z^ Z rn . < 1 defines a domain of validity 
for the coupled system of Dyson-Schwinger equations (f49|) and (|50| ) and the solution 
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of this system of equations gives the propagators that are to be used in the BS 
equation (f0|). Therefore the condition < Zj, Z m . < 1 gives an upper bound for the 
coupling constant g above which the quanta of the field $2 n o longer correspond to 
physical particles and above which the ladder-approximation to the BS equation is 
not applicable. 



4 Bound States in QED 



It is certainly interesting to compare the results obtained so far in the scalar model 
with an exactly renormalizable theory. QED whose Dyson-Schwinger and BS equa- 
tions have been studied intensiveley provides hereby a good testing ground for 
interpreting our results. We first calculate the spectrum of the ladder approximation 
to the BS equation for positronium, and we will compare to the results that have 
been obtained for the scalar theory. As in the last section we will derive consistent 
Dyson-Schwinger and BS equations, taking one approximation for the CJT-action 
as starting point. A simplified version of the coupled system of Dyson Schwinger 
equations will be solved. 

First, we will discuss the BS equation in ladder approximation for positronium. 
However, we will also present results for constituents of different mass and for a 
massive photon when commenting on the problem of complex eigenvalues. The ladder 
BS equation is given by 

T (q, P) = -ie 2 [ (q - k) 1 ,S 1 (k+) T (k, P) S 2 (*_) lv (53) 

J (2vr) 

k + = k + aP 

k- = k — (l — a)P 

where V is the BS vertex function which is related to the BS amplitude as given in 



fllOD . The propagators of the constituents are approximated by the corresponding 



bare propagators 

Si(p)= J + T\- > i = U. (54) 

D^ v is the bare propagator of the photon, and since we are working in the ladder 
approximation we used a bare fermion-photon vertex, i.e. 7^ only. 

The general expression for the BS vertex function T describing a pseudoscalar 
bound state was derived in [E3|. However, in the Feynman gauge and considering 
only the ladder approximation one can derive that the tensor part of the vertex 
function vanishes. Taking this into account we arrive at 

T (q, P) = 75 (ri (q, P) + {P-q)i Y 2 (q, P) + P T 3 (q, P)) (55) 

19 



200- 



150- 



100 



50- 



• • • 

• • • 

• • • 



iiiiili 



— i — 1 — i — 1 — i — 1 — i — 

20 40 60 80 



"0.807 



300 



number of mesh points for the momentum integration 



Figure 10: Solutions of the BS equation (j54|) for constituents of equal mass in the 
limiting case M —>■ 0. The solutions are given for < a < 200 and for various 
numbers of meshpoints. 



Substituting this ansatz into the BS equation (|54]) and taking appropriate combina- 
tion of traces and projections one arrives at a coupled system of integral equations 
for the scalar functions r 1; T 2 and T 3 . The details of the numerical method are given 
in appendix [D], and we proceed immediately to the discussion of the results. 

First, we consider the spectrum for constituents of equal mass in the limit of 
massless bound states. Goldstein [2B[ was the first who considered this limiting 



case. He found massless solutions for all values of the coupling constant. One can 
demonstrate that the BS equation ([54]) in this limit reduces to an equation for Ti 
which is not coupled to r 2 and r%. Therefore one may impose the "normalizibility 
constraint" 

J d 4 k\T 1 (k)\ < oo. (56) 

Taking this condition into account Goldstein still found a continuum of solutions 
which, however, starts at a = e 2 /4n = tt/4. The results of the numerical solution 
are displayed in fig. ID for an increasing number of mesh points for the numerical 
momentum integration. Hereby every point represents a solution of the BS equation 
for M = (within numerical accuracy), and the corresponding value of the fine 
structure constant a is represented. As the number of mesh points is increased one 
finds an ever increasing density of solutions in the intervall < a < 200. This 
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Figure 11: Solutions of the BS equation fl54|) for constituents of equal mass. The 
parameters are defined as a = e 2 /4n and t] = M/(m\ + m 2 ). The lines represent 
discrete solutions and the shaded area stands for a continuum of solutions. We took 
40 meshpoints for the momentum-integration which lead to a lower bound of a ~ 0.59 
for the continuum. 



is exactly the way a continuum of solutions is expected to show up in a numerical 
treatment. However, the lower bound of this continuum depends on the number of 
mesh points for the momentum integration. Taking 40 meshpoints the continuum 
starts at a ~ 0.59 and taking as much as 300 meshpoints we found the lower bound 
to stabilize at a ~ 0.807. This compares reasonably well to a = 7r/4 0.785 which 
was found by Goldstein. Note that this continuum of solutions is not related to the 
physical continuum of the scattering states. 

We turn to the case of rather weakly bound states of the electron and the positron 



(weakly bound equal mass case). The results are shown in fig. |TT|. The lines represent 
discrete solutions which start at a ~ as well as at greater values of a and the shaded 
area marks the beginning of a continuum of solutions. Since we took 40 meshpoints 
for the standard calculations the continuum starts at a ~ 0.59. 

For < i] < 0.94 the continuum is the only feature of the spectrum of solutions 
of the BS equation ([54]), i.e. the continuum of solutions (which was known to appear 
for M — > 0) is present for all values of r], and it is so with a constant lower bound. 



21 



The BS equation for scalar constituents and scalar exchange particle gave complex 
solutions at(j]) for constituents of different mass and for a massive exchange particle. 



We solved the BS equation (54) for various mass ratios m.i/m.2 7^ 1 and m 7 /m2 7^ 
and the results agree qualitatively with those that have been obtained for the scalar 
theory. To mention two of the more pronounced features: a massive photon destroys 
the clustering of solutions at a = and for M — > m\ + mi. As well as for the scalar 
theory we found complex solutions 01(77) which, however, were found only above the 
lower bound of the continuum. 

Let us finally take up the discussion of the (in-)consistency of the approximations 



for the propagators and for the kernel that enter the BS equation (|54|) . As for the 
scalar theory we approximate the loop series T 2 by 



2 



J d A xd 4 ytr(s(x,y)^S(y,x)YD ll A^y)) (57) 

where S and D^ u are the propagators of the electron and the photon, respectively. 
These propagators are not the corresponding bare propagators but the solution of 
the coupled system of Dyson- Schwinger equations: 

5T 5T 

«r«v= a (58) 

Using the approximation (^) for T 2 one may derive the explicit Dyson-Schwinger 
equations that correspond to (|58|). This leads to a system of equations which is shown 



in fig. 12 and which compares to the system of equations for the scalar theory fig. 



If the appoximation ([)?]) is used as input for the general form of the BS equation 
dV dV rg> (x, y- x\ y') r' k ' (x f , y') = (59) 



one obtains the ladder approximation to the BS equation which is given in 
However, it has to be emphasized that the propagators S and D are not the bare 
ones but the solution of flSlf). 

For the scalar theory we showed that the physical condition < Z iy Z m . < 1 
gives a domain of validity for the system of Dyson-Schwinger equations as well as for 
the ladder approximation to the BS equation. In order to obtain a first estimate we 
solved the system of equations in quenched approximation thus effectively neglecting 
the equation for the photon. The equation for the propagator of the electron is then 
given by 

S- 1 (p) =Z 2 (tf- m ) - iZ,e 2 T -^fS (k) YD» U (k - p) (60) 

J (2tt) 
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Figure 12: Graphical representation of the coupled system of Dyson-Schwinger 
equations (|58|). Thick and thin lines represent the dressed and the bare propagators, 
respectively. 

where the bare mass mo as well as the renormalization constants Z\ and Z 2 depend 
on the renormalization scale \i and on the cutoff A. As one can read off from ( j60|) the 
approximation (|57|) yields the bare fermion-photon vertex together with an obviously 
nonperturbative propagator for the fermion. This implies that the Ward-Takahashi 
identity cannot be satisfied, i.e. Z x = Z 2 cannot be assumed. Since the use of the 
bare vertex enforces Z\ — 1 we determined only Z 2 and m/m according to the 
renormalization condition 

S~ 1 {p)\ p 2 =0 =lt>-m (61) 

and fixed Z\ to 1. For simplicity we employed the normalization at the soft point 
p 2 = instead of the on-shell renormalization. Defining the scalar functions A and 
B through the relation 

s{p) = pWW^bW) (62) 



one obtains from (|60|) a coupled system of equations for A and B after taking appro- 
priate combination of traces. In terms of these scalar functions the renormalization 
conditions read A(0) = 1 and B(0) = m. 

The details of the numerical solution are given in appendix [E| and therefore we 
focus on the results. For the physical coupling constant a = 1/137 the scalar functions 
A and B are nearly equal to their perturbative limits, i.e. A ~ 1 and B w m, as could 
be anticipated. 

For greater values of a an increasing number of iterations was needed in order to 
obtain a self-consistent solution. We found no solution for a > 0.74 which compares 
quite well with the critical value a = 7r/4 w 0.785 and with the lower bound of the 



continuum, which was estimated to be at a ~ 0.807. Fig. 13 displays the mass-ratio 
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Figure 13: The mass-ratio m/rriQ versus the coupling constant a. The divergence 
of this quantity signals chiral symmetry breaking. 

m/rriQ. The divergence of this quantity signals dynamical generation of a mass for 
vanishing bare mass, i.e. spontaneous breaking of chiral symmetry. To clarify the 
significance of this result in comparison to the scalar model we plot in fig. |14] the 
inverse mass-ratio mo/m and the renormalization constant Z 2 versus the coupling 
constant a. Both Z 2 and mo/m change sign at the same value of a, a ~ 0.75. A 
commonly accepted interpretation of this is that dynamical chiral symmetry breaking 
occurs at the critical value of the coupling constant, and that above this value of the 
renormalized fine structure constant QED is not defined [0, [28[ . 



5 Conclusions 

The central aim of this investigation was a clarification how (or more precisely, 
whether) the physical spectrum of excited states may be infered from the fully rela- 
tivistic Bethe-Salpeter equation. In a very strict sense the answer is negative: Taking 
into account the need for renormalization in the underlying theory the applicability 
of the ladder Bethe-Salpeter equation turns out to be very restricted. The com- 
parison to (quenched) QED demonstrates this clearly. It is asserted that the value 
of renormalized fine structure constant has to be lower than some critical value in 



order to allow for reasonably behaved quantum field theory, see e.g. |28] for a cor- 



responding result in a Dyson-Schwinger approach and |29| for a lattice calculation. 
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Figure 14: The renormalization constant Z 2 and the mass-ratio m /m versus the 
coupling constant a. 

Finding the similar behaviour that the field renormalization constants becomes neg- 
ative for large couplings one may speculate that our (super-renormalizable) scalar 
model is also well-defined only for renormalized couplings below the critical one sig- 
naled by the zero of the renormalization constant. Below these critical values of 
the coupling constant the Dyson-Schwinger equations for the propagators as well as 
the Bethe-Salpeter equation in rainbow-ladder approximation were shown to give 
physically acceptable solutions only. In addition, the use of bare propagators in the 
Bethe-Salpeter equation introduces only a small additional error. 

For large couplings the homogeneous ladder Bethe-Salpeter equation possesses 
abnormal solutions being excitations in the relative time coordinate. Half of these 
abnormal states are odd w.r.t. this relative time, and they lead to complex eigenvalues 
via the crossing with even states. In addition, in QED there is the problem of a 
continuum of solutions (Goldstein problem). The investigation reported here provides 
evidence that all these crazy things happen beyond the applicability region of the 
underlying renormalized quantum field theory. Thus, before studying a spectrum 
using the Bethe-Salpeter equation one should ensure oneself that the renormalization 
of one-particle properties can be performed with a physically acceptable result. 
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Appendix 



A Numerical method for the solution of the scalar 
Bethe— Salpeter equation 

The Wick-rotated homogeneous BS equation for scalar constituents and scalar ex- 
change particle is given by 

( P\ = 2 1 1 f d 4 p' X (p', P) f . 

X & ^ 9 v\ + m\ pi + m 2 J (2vr) 4 (p - pf + ^ ' [M) 

where the euclidean momenta pi are defined according to 

px-.— aP+p, p% :— (1 — a)P — p. (64) 



The symmetries of this equation have been discussed in section |0| and are sum- 
marised in table [|. For a massless exchange particle (|63|) is invariant under 0(4) 
transformations and the BS amplitudes are proportional to a spherical harmonic of 
SO (A), denoted by Zki m . In the case of 0(3) angular momentum I = the Zkim 
are proportional to Gegenbauer polynomials Cl of degree 1, see |3(J for their def- 
inition. For a massive exchange particle the 0(4)-symmetry is only approximate. 
Nevertheless, it turns out that the expansion 

CO 

X (p, P) = Y,C 1 m (cos (9)) xm (INI , ||P||) 9 := < (p, P) (65) 

m=0 

still converges very fast. For further usage we now define the dimensionless "mo- 
menta" 



x := i^l x' := , (66) 
m 2 m 2 

the dimensionless mass-parameters 

f:=I5l, K:= A „:=_" (67) 
m 2 m 2 mi + m 2 

and the dimensionless coupling constant 
167r z m 2 

Furthermore the abbreviations 

K 2 + X 2 + X' 2 



t , B := z - Vz 2 - 1 (69) 
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will be used. The integration over the euclidean momentum p' is done in spherical 
coordinates: 

d A p' = dp'p' 3 de'sm 2 6' d(p' simp' & := < (p', P) .(70) 

Jo Jo Jo Jo 

The integrations over 9',(p',ip' can be done analytically and the integral over the 

absolute value of p' has been mapped to the interval [—1, +1] according to 

a/ = tan(-(l + 2)) !°° dx' ^- [ + (71) 

l 4 l "J, 47_! cos 2 (f(l + ^)) 1 J 

for the numerical ^-integration 

+1 Nl 

dzf(z) -> £>,(;)/(*,). (72) 
1 i=i 

w z (j) are hereby the weights of the Gauss- Legendre integration. After projection 
onto the expansion coefficients Xk we arrive at an eigenvalue problem for the coupling 
constant A 

<\>k fa, 1]) = A ^2 ^ ' 0" ^) ' ( 73 ) 

n=0 3=1 

where the auxiliary functions 0^ := i k x 2 Xk{p, P) have been used. The indices n and 
z run over the degree of the Gegenbauer polynomial and over the Ni mesh points for 
the momentum integration, respectively. The definition of the matrix n j is given 
by 



+ 1 

5 (1 - 7] 2 ) + x 2 1 - rf + x 2 



fel,ry ' 4 <V*i ^ + y 2 cos 2 (f (1 + Z jj) n + 1 ' 1 J 



" V y V ' (75) 

' 2 ()//./•, " 2//./-, 

(-1)" ^ (_i)^+ fc +"4"- fc l+ 1 + ^ +fc+3 

+ T+Af —J' (77) 

(7 := ^(|n-fc| + A; + n), (78) 

and (79) 



Ai := Vi-y/yf + l. (80) 

Different LAPACK-routines have been used to solve the eigenvalue problem (|73|) and 
to cross-check the results. 
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B The Bethe— Salpeter amplitudes for euclidean 
momenta an for constituents of equal mass 

The BS amplitudes describing a bound state of scalar constituents of equal mass are 
defined by 

X {x 1 ,x 2 ,P) = (0\T{^{x 1 )^{x 2 )}\P) (81) 
x(x h x 2 ,P) = (0\T{&( Xl )&{x 2 )}\P) (82) 
= (0\T{^(x 1 )^(x 2 )}\P)* 

where T is the usual time ordering operator and T denotes time ordering in the 
reverse order. Now we use the definitions 

X = axi + (1 — a)x 2 (83) 
x = Xi — x2 

to separate the center of mass motion 

X (x 1 ,x 2] P)=e- lP - x 4>(x,P) (84) 

defining thereby <f), which denotes the relative part of the BS amplitude. Taking this 
separation of the center of mass and writing the time ordering explicitely we obtain 

X ( Xl ,x 2 , P) = e~ lP - x (0 (a; ) g (x, P) + (-x°) g (-x, P) ) (85) 
x(xi,x 2 ,P) = e* px (Q(x°)g(-x,P) + Q (-x°) g(x,P))* . (86) 

Since we assumed constituents of equal mass we may take the amplitudes to have a 
definite x -parity e (0) which is defined by 

0(-x°,x) =e(0) 0(x°,x) . (87) 

This can be used to rewrite eqs. (|85| ) ,(|86| ) in the following way 

( Xl ,x 2 , P) = (0 (x°) g (x, P) + g (-x, P) ) , (88) 

0(x u x 2 ,P) = e(0) (e(x°)g(-x°,x,P) + & (-x°) g (x° , -x, P))* . (89) 

Using the relation 

(z) = — — f dke~ ikz , 1 - (90) 
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we find the following Fourier-transforms 

/ <* 4 --e(z») ff (*,P) - ~ J ^f^r + < (»> 
/A e -"«e(-x°)[ ff (x°,- x ,p)]* = ~ [dq' ^f'V;^ 

J /J 2 717 J g u +p u + zO+ 

which may be used to rebuild the Fourier transforms of the relative part <fi of the BS 
amplitudes. After Wick-rotation p — > ip^ we find the relation 

^(p,P)=e(0)[j(p,P)]* (92) 
for the relative parts of the BS amplitudes and for euclidean relative momenta p. 

C Numerical method for the solution of the Dyson- 
Schwinger equations of the scalar theory 

The coupled system of DS equations for the propagators of the scalar fields $j , i — 
1, 2, 3, is given by 

—l I G 3 (p-q)G i (q) , i = l,2, (93) 
(2vr) 

S 3 (p 2 ) = A 3 + g 2 i I £ Z 4 G, (p - g) (g) (94) 

7 (2tt) . =1 

where the abbreviation 

A, := (1 - Zi) p 2 + {ZiZ m . - 1) m\ , % = 1, 2, 3 , (95) 

and the definition £j(p 2 ) := p 2 — m 2 — G^ 1 (p) have been used. For notational conve- 
nience we introduce 

—^-G t (p-q)G k (q), i, k =1,2, 3, (96) 

which allows to rewrite the equations fl9~3"D,(|94]) in the form 

Si(p 2 ) = (l-^)p 2 + (^Z mi -l)m 2 + Z 4 2^P 3i (p 2 ), * = 1,2, (97) 

2 

S 3 (p 2 ) = (l-^^+^ma-^^+Zi^Z^Cp 2 ). (98) 



i=l 
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We use the on-shell renormalization condition £,; (m 2 ) = to eliminate Z m . and 



obtain 

E,(p 2 ) = (l-Z l )(p 2 -m^ + Z l 2g 2 t(F 3t (p 2 )-F 3l (m 2 i )y i = l,2, (99) 

2 

£ 3 (p 2 ) = (l-Z 3 )(p 2 - m 2 )+g 2 iJ2Zi(mP 2 )-FiM)). (100) 



i=l 



The contain the self-energies £j ((p — g) 2 ) and £& (q 2 ). In the following we use 
an angle approximation, i.e. we substitute max(g 2 ,p 2 ) for (p — q) 2 if (and only if) 
(p — q) 2 appears as the argument of £j , i — 1, 2, 3. After Wick-rotation 

9 = (0o, q) -»• (*?4, q) =: (101) 
we change to spherical coordinates 

rAuY fit /*27T 

d 4 q^ / (i00 3 / d6 sin 2 6 / dcpsincp / dip (102) 
Jo Jo Jo Jo 

where we temporarily introduced the cutoff A uv . Because of the angle approximation 
the integrations over 9, if and ip can be done analytically which leads to 



ik 



2 v^c 2 + 4p 2 g 2 - c 
rf 9 2 , ^ ? =^ > ( 103 ) 



32ttV 7 <f + m 2 + £ fc (-<f, 
where the abbreviation 

c := -p 2 + q 2 + m 2 + E< (- max(p 2 , 2 )) (104) 

has been used for convenience. The system of equations has been solved iteratively 
and in each step of the iteration the field renormalization constants Zj are calculated 
according to 



ap z 



0. (105) 



The solutions , i — 1, 2, 3 , of the system of equations (|93| ) and ( |94"|) enter the 
BS equation via the propagators of the constituents and the exchange particle, see 
equation (f44"|) and its derivation. More explicitely: 



^ (p, P) = agH Gl (p0 G 2 (p 2 ) / , 2 f 2 } q ;l\ ( ~?K (106) 

J (2tt) (p-g) +m| + £ 3 ((p-g) ) 

where 

Pt:=p-aP, p 2 :=p+ (l-a)P (107) 
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and 

G < W = i?T^W)' i = 1 ' 2 ' (108) 

This equation is no longer an eigenvalue problem for the square of the coupling 
constant g 2 . However, one may introduce a formal eigenvalue a and regard g and 
M as parameters. A solution, i.e. a pair (g, M) is now signaled by the eigenvalue 
a = 1. Since we use the angle approximation it is sufficient to apply the substitution 
mf — > mf + Sj(p^ ax ) to the eigenvalue problem as given in (|73|) and the following 
equations in order to obtain the corresponding eigenvalue problem with self-energies 
taken into account. 



D Numerical method for the solution of the Bethe— 
Salpeter equation for QED 

The BS equation for positronium is given by 

r (q, P) = -ii / d A k (q - k) ^G, (k + ) T (k, P) G 2 (&_) 7, • (109) 
(2tt) J 

The propagators of the constituents and the corresponding momenta are defined 
according to 

n(\ / + m i k + :=k + aiP , , 

Gi{p) = ^ 2 , , , , , ai + a 2 = l , (110) 

p 2 _ m 2 + u k_:=k- a 2 P 

where we allowed for different masses of the electron and the positron. Since all the 
calculations have been done in the Feynman gauge we used 

JV(P) = -TTT-, (I") 
p z + ie 

T (q, P) = 75 (r\ (q, P) + {P- q)^T 2 (q, P) +/T 3 (q, P)) , (112) 

for the propagator of the photon and for the vertex function. (For discussion and 
references see section |j.) As for the scalar theory we use Gegenbauer expansions for 
the scalar functions 

oo 

r*(g ) P)=^C?(cose)[r fc ] i (||g||,||P||) , 9:=<(q,P) , (113) 

i=0 

as well as for the propagators of the constituents and the photon. Applying appropri- 
ate combinations of traces and projections one obtains a coupled system of equations 
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for the coefficients [r^ of the Gegenbauer expansion ( |lf3| ). After Wick-rotation 
V — (Po, p) - > {ip4, p) = : P this system of equations reads 

d 3 yP t Y,C!,(cos6 k )S 3 S 2 S 1 (111) 

i'=0 



_j? 

(27r) 3 axa 2 + I) 2 xif 



(A n [Tx],, + a 12 [r 2 ], + a u [r 3 ]. 



i'J ' 



-1 



(27r) ai«2 (5 + 1) ^ 2 



r 00 1 

/ dPyP i ^2cl,(cose k )S 3 S 2 S 1 -=- (115) 



( [i (8 + 1) (rjx cos9) A 21 - (6 + if V 2 A 31 ] [T 1 ] l , 
+ [i (5 + 1) (rjx cos 6) A 22 - (5 + l) 2 r/ 2 A 32 ] [T^, 
+ [i (5 + 1) (77a; cos 0) A 23 - (5 + l) 2 r/ 2 A 33 ] [r 3 ],,) 



/ dPyPi^CKcoaejSaS&y (116) 



(2^)^ ai« 2 (5 + l) 2 XT] 2 „ 

( [-(5 + lfr 1 2 A u ~x 2 A 21 ] pi],, 
+ [- (8 + lfrj 2 A 32 -x 2 A 22 \ \T 2 ] V 
+ [-(Hl) 2 i| 2 4-x 2 4] [T 3 ] lf ) 
In addition to the dimensionless mass-parameters 



M 



V 



mi 



mi + Tn>2 m 2 
and the dimensionless "momenta" 

Ikll ||*|| 

m 2 m 2 
several abbreviations have been used: 
2 



7T 



dO sm 2 (8)Cl(cos6) 



PiCl(cos6) = 5 t 



d 3 y 



dy / d9 k sin 9 k / dip k sin ip k 
Jo Jo 



Si 



S 2 



91 



11=0 



n 2 =0 

CX) 



#2 



ni+1 



71 3 =0 



#3 



n 3 +l 



n 2 +l 



(117) 
(118) 

(119) 
(120) 

(121) 
(122) 
(123) 
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. y 2 -a 2 (£ + l)V + l . y 2 -a 2 (5 + l) 2 rf + <5 2 

2a 2 {£> + 1) yr] 2an (d + 1) yr] 

= xf_+y^ , , 

ys 2xy ' V ; 



/ 2 = — i (5 + l) 3 (r]xcos8) [(r]xcos9) 2 + r^x 2 ] , (126) 
f 3 = (5 + I) 2 {-x 2 (5 + l) 2 rf + (r 1 x cos 6f) , 

A n = 4( 2 / 2 -z(a 1 -a 2 )(5 + l)(wcos4) + aia 2 (5+l) 2 r/ 2 + 5), (127) 
A 12 = 4i(5+l)( W cos^)((l-5)y 2 -i(a 1 + a 2 5)(5 + l)( W cos^)), 
A 13 = A (i (5 2 - l) (r]y cos 9 k )- (a, + a 2 5) (5 + l) 2 r] 2 ), 

A 21 = -2i{5 2 -\) (r 1 ycos9 k ) + 2(a 1 + a 2 5)(5 + l) 2 rf, 
A 22 = i (6 + 1) (r]y cos 9 k ) (2r (-y 2 + a,a 2 (5 + l) 2 r] 2 -5) (5 + 1) (r]y cos 9 k ) 

+2 {a 2 - ai )y 2 (5+l)V 
-Aiaia 2 (5 + l) 3 rf (rjy cos 9 k )) , 
A 23 = -4(5+l) 2 (r ] ycos9 k ) 2 + 2(5 + l) 2 y 2 r ] 2 

+2i («i - a 2 ) (5 + l) 3 rf (rjy cos 9 k ) 
-2a x a 2 (b + \frf-2b (5 + l)rf, 

A 3 i = 2 (5 — 1) (xy coscj) + 2i (a± + a 2 5) (xr/ cos 9) (5 + 1) , 
A 32 = i (5 + 1) (rjy cos 9 k ) (-2 (-y 2 + a x a 2 (5 + I) 2 rj 2 - 5) (xy cos uj) 

+2i [a 2 — af) (5 + 1) y 2 {xr] cos 9) 
+ 4ai<y 2 (5 + l) 2 (rjy cos 9 k ) (xrj cos 9)) , 
A 33 = —Ai (5 + 1) (-qy cos 9 k ) (xy cos uj) + 2i (5 + 1) y 2 (xr] cos 9) 

—2 («x — a 2 ) (5 + l) 2 rf (xy cos uj) 
—2%a\a 2 (5 + l) 3 rf (xrj cos 9) 
-2i5 (8 + 1) (xrj cos 9) , 

9 := < (q, P) , 9 k :=< (k, P) , uj := < (q, k) . (128) 



As for the scalar BS equation we used the transformation (71]) to map the integra- 
tion over the absolute value \\k\\ to the range [— 1,+1]. The integrations over \\k\\ 
(corresponding to z after the transf. (|7T|) ) and over the angles 9 k ,<p k ,9 have been 
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performed with the Gauss-Legendre quadrature. Finally one obtains an eigenvalue 
problem for the square of the coupling constant which can formally written as 



9 





B\2 


Bn \ 






Bn 


B22 


B23 


II 


r 2 


B31 


B32 


B33 J 




^ r 3 



(129) 



The B^ are matrices of dimension Ni (N 2 + 1) x Ni (N 2 + 1) where Ni is the number 
of mesh points that have been used for the momentum integration over z and A^ is 
the maximal degree of Gegenbauer polynomials that have been taken into account. 

E Details of the numerical solution of the Dyson- 
Schwinger equation for the electron propagator 

Using the definition 

S(p)= ^A^)-B(jP) (130) 
of the scalar functioins A and B the DS equation (|60| ) translates into 

A{f)f + B{f) = Z 2 (j/-m ) (131) 

" A d A k (-1A(k 2 )Y + AB{k 2 



(2tt) 4 V ^4 2 (k 2 ) k 2 - B 2 (k 2 ) J (p - kf 



Applying |tr ( ) and |ir (7^ ) leads to separate but coupled equation for the scalar 
functions A and B. After Wick-rotation p — (po, p) — > (ip4, p) ='■ P the coupled 
system of equations for A and B reads 

+ wf* s wlp) (132) 

+ ?4r"at- - 4( - p) 



8tt 2 J P 2 A 2 (-k 2 ) k 2 + B 2 (-k 2 ) ' 

B(-f)=z 2mo + d »3__^_i_ (133) 



2tt 2 J P 2 A 2 (-k 2 ) k 2 + B 2 (-k 2 ) 
where now all the momenta are euclidean. 
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In order to derive eqs. 



(p-k)' 

p ■ k 
(p-k)' 



2vr k 



7T~ 



and ( |133| ) we used the relations 

e(* 2 -p 2 )pH-e(p ! -P)i 
e(P-p 2 )i + e(p 2 -P) * 2 



(134) 
(135) 



where 



/pn pit 1*211 

dVL A := / d9k sin 2 9^ / dcff. sin^ / 1 
Jo Jo Jo 



dip. 



k ■ 



(136) 



The wave function renormalization constant Z 2 and the mass ratio m/m have 
been fixed according to 



A(0) 



5(0) = m. 



(137) 



and the vertex renormalization Z\ has been set equal to 1, for a discussion of this 
point see section |j. 

We took an exponential distribution of mesh points ki and used an extended 
formula of order 0(1/N 3 ) to perform the momentum integrations, see [plj . 
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